Formal matched asymptotics for 
degenerate Ricci flow neckpinches 



Sigurd B. Angenent, James Isenberg, and Dan Knopf 

Abstract. Gu and Zhu 1161 have shown that Type-II Ricci flow singulari- 
ties develop from nongeneric rotationally symmetric Riemannian metrics on 
S™+ 1 (n > 2). In this paper, we describe and provide plausibility arguments 
for a detailed asymptotic profile and rate of curvature blow-up that we predict 
such solutions exhibit. 



Let (§™ +1 , g(t) : < t < T) be a rotationally symmetric solution of Ricci flow. 
Gu and Zhu prove that Type-II singularities can occur for special nongeneric initial 
data of this type |16j . Garfinkle and one of the authors provide numerical simula- 
tions of the formation of such singularities at either pole [131 114) . However, almost 
nothing is known about the asymptotics of such singularity formation, except for 
complete noncompact solutions on M 2 , where Ricci flow coincides with logarith- 
mic fast diffusion, u t — A log u. (The asymptotics of logarithmic fast diffusion, 
which are unrelated to the results in this paper, were derived by King |18) and 
subsequently proved in K 2 by Daskalopoulos and Sesum |10|.) 
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We say a sequence {(pi,ti)}^ of points and times in a Ricci flow solution is 
a blow-up sequence at time T if t , /* T and |Rm(a;j,tj)| — > oo as i — > oo. We say 
(M"" +1 , g(£)) develops a neckpinch singularity at T < oo if there is some blow-up 
sequence at T whose corresponding sequence of parabolic dilations has as a pointed 
limit the self-similar Ricci soliton on the cylinder R x §™. We call a neckpinch 
nondegenerate if every complete pointed limit formed from a blowup sequence at T 
is a solution on the cylinder; we call it degenerate if a complete smooth limit of some 
blow-up sequence has another topology. Rotationally symmetric nondegenerate 
neckpinches have been studied by Simon [21| and two of the authors [H [3], who 
proved they are Type-I ("rapidly forming") singularities in which the curvature 
blows up at the natural parabolic rate, (T — t) sup pgM „+i |Rm(p, t)\ < oo. On 
the other hand, Type-II ("slowly forming") singularities have the property that 
(T — t) sup peM n+i |Rm(p, t)\ = oo. The compact Type-II singularities proved to 
exist by Gu and Zhu |16j are degenerate neckpinches. 

This paper is the first of two in which we study the formation of degenerate 
Ricci flow neckpinch singularities. In the present work, we assume that a degenerate 
neckpinch singularity occurs at a pole at time T < oo, and we derive formal matched 
asymptotics for the solution as it approaches the singularity. This procedure pro- 
vides evidence for a conjectural picture of the behavior of some (not necessarily all) 
solutions that develop rotationally symmetric degenerate neckpinch singularities. 
In particular, it predicts precise rates of Type-II curvature blow-up^] In forthcom- 
ing work, we will provide rigorous proof that there exist solutions exhibiting the 
asymptotic behavior formally described here. 

Given any rotationally symmetric family g(t) of metrics on S n+1 , one may 
remove the poles P± and write the metrics on §™ +1 \{-P±} ~ (—1,1) x §" in the 
form g(t) — (ds) 2 +if) 2 (s,t) <7 can , where s(x,t) denotes g(i)-arclength to x £ [— 1, 1] 
from a fixed point xq £ (— 1, 1). (Here, g can is the canonical unit sphere metric on 

see Section [2] for a detailed discussion of these coordinates.) Thus the function 
ip(s(x,t),t) completely characterizes a given solution. Our Basic Assumption in 
this paper, explained in detail in Section[2j is that the initial metric g(s, 0), hence 
ij)(s, 0), satisfies certain curvature restrictions which ensure that the geometries we 
consider are sufficiently close to those studied in [13 , !|14] and [16] , and therefore 
are likely to develop neckpinches. It also allows us to employ certain prior results 
of two of the authors [2], which are useful for the arguments made here. (Compare 

P30 

In the language of Section [2] let s(t) denote the location of the local maximum 
("bump") of i/j closest to the right pole x = 1. Note that s(t) may be only an 
upper-semicontinuous function of time, because a bump and an adjacent neck can 
join and annihilate each other. Lemma 7.1 of [2] proves that the solution exists 
until tp becomes zero somewhere other than at the poles. Since there is always at 
least one positive local maximum of ip, the quantity s(t) is defined for as long as 
the solution exists. Lemma 7.2 of [2] proves that if lim^Ti tp(s(t), t) > 0, then no 
singularity occurs at the right pole. Hence, if a degenerate neckpinch does happen 
at the right pole, it must be that Hmty<Tip{s(t),t) = 0. This can happen only if (i) 
the radius tp vanishes on an open set, or (ii) the bump marked by s moves to the 
right pole. Note that these alternatives are not mutually exclusive. In either case, 



For a comprehensive statement of our results, see Section 7 
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it follows that there are points s(t) such that 

(1.1) lim s(t) = s(l,T) and lira tp s (s(t),t) = 0. 



Observe that (1.1) is incompatible with the boundary condition (2.3) necessary 
for regularity at the pole, which forces ip = (s(l) — s) + o(s(l) — s) as s /• s(l). 
Therefore, a consequence of our Basic Assumption is that if a degenerate singularity 
does develop, then this expansion for ip cannot hold uniformly in s as t f~ T. 
Instead, the solution must exhibit different qualitative behaviors in a sequence of 
time-dependent spatial regions. 

Starting with that fact, we construct in this paper a conjectural model for 
rotationally symmetric Ricci flow solutions that develop a degenerate neckpinch, 
and we check its consistency. We do this systematically by studying approximate 
asymptotic expansions to the solution in four connected regions, which we call 
(moving in from the pole) the tip, parabolic, intermediate, and outer regions. By 
matching these expansions at the intersections of the sequential regions, we produce 
our model. This is the essence of the formal matched asymptotics process familiar 
to applied mathematicians. The process involves first formulating an Ansatz, which 
consists of a series of assumptions (listed below) pertaining to the geometry and 
its evolution equations in the four regions. Then, matching across the boundaries 
of the regions, one constructs formal (approximate) solutions. Finally, one checks 
that these formal solutions remain consistent with the Ansatz, and argues that the 
approximations remain sufficiently accurate. 

For clarity of exposition, we establish notation in Section [2] and then begin our 
study working in the parabolic region rather than the tip. We treat the tip (the 
most critical region) last. Finally, in Section [7j we provide a detailed summary of 
our results and the conjectural picture they provide. 

2. Basic equations 

We begin by recalling some basic identities for SO(n + l)-invariant Ricci flow 
solutions. 

To avoid working in multiple patches, it is convenient to puncture the sphere 
S n+1 at the poles P ± and identify S n+1 \{P±} with (-1, 1) x S". If we let x denote 
the coordinate on the interval (—1, 1) and let <7 can denote the canonical unit sphere 
metric, then an arbitrary family g(t) of smooth SO(n+ l)-invariant metrics on § n+1 
may be written in geodesic polar coordinates as 

(2.1) g(t) = <p 2 (x, t) (dx) 2 + xjj 2 {x, t) 3can . 

Denoting the distance from the equator {0} x S™ by 

s{x,t) = f <p(t,t)d£, 



allows one to write (2.1 1 in the more geometrically natural form 
(2.2) g = {dsf +^ 2 {s{x,t),t)g C3 . n , 

where "d" is the differential with respect to the space variables but not the time 
variables. We also agree to write 

d_ _ 1 d 

ds f{x, t) dx 
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With this convention, d/ds and d/dt do not commute; instead one gets the usual 
commutator 

d d~\ ft d 

dt' ds\ ip ds 

Smoothness at the poles requires that ip satisfy the boundary conditions 
(2.3) lU a= ±i=Tl. 



The metric (2.2 1 has two distinguished sectional curvatures: the curvature 

ip 2 



( 2 - 4 ) L = -,,2 



of a plane tangent to {s} x S™, and the curvature 

(2.5) K = -^ 

ip 

of an orthogonal plane. 
Under Ricci flow, 

h = - 2Rc ' 

the quantities ip and ip evolve by the degenerate parabolic system 



ft = n 



cpip ip 2 ip 



ipt = —2 5- + (n- lj-jy — . 

ip z ip 3 ip z ip ip 

The degeneracy is due to invariance of the system under the infinite-dimensional 
diffeomorphism group. By writing the spatial derivatives in terms of d/ds, one 
can simplify the appearance of these equations, effectively by fixing arclength as a 
gauge. The quantities ip and ip in (2.1) then evolve by 

(2.6a) tpt = n^-<p, 

ip 

(2.6b) ^ = ^ s _( n _i)izA 5 

ip 

respectively. 

In [2] , nondegenerate neckpinch singularity formation is established for an open 
set of initial data of the form ( 2.2 1 on S" +1 (n > 2) satisfying the following assump- 
tions 

(1) The sectional curvature L of planes tangent to each sphere {s} x S" is 
positive. 

(2) The Ricci curvature Rc = nK(ds) 2 + [K + (n — l)L]ip 2 g can is positive on 
each polar cap. 

(3) The scalar curvature R = 2nK + n(n — 1)L is positive everywhere. 

(4) The metric has at least one neck and is "sufficiently pinched" in the sense 
that the value of ip at the smallest neck is sufficiently small relative to its 
value at either adjacent bump. 

In [2, precise asymptotics are derived under the additional assumption: 



2 We call local minima of ip "necks" and local maxima "bumps" . The "polar caps" are the 
regions on either side of the outermost bumps. 
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(5) The metric is reflection symmetric, tp( s ) = V-K - s )i an d the smallest neck 
is at x = 0. 

For our analysis here of degenerate neckpinch solutions, we start from initial 
geometric data that satisfy some but not all of these assumptions. In particular, 
we drop the condition on the "tightness" of the initial neck pinching, and we also 
do not require reflection symmetry^] We now state our underlying assumptions on 
the initial geometry: 

Basic Assumption. The solutions of Ricci flow we consider in this paper are 
SO(n + l)-invariant. They satisfy conditions (l)-(3) above, and they have initial 
data with at least one neck. Furthermore, we assume that a singularity occurs at 
the right pole (x = +1) at some time T < oo. 

3. The parabolic region 

A consequence of our Basic Assumption is that R > at t — 0. By standard 
arguments, this implies that the solution becomes singular at some time T < oo. 

We assume that a neckpinch occurs at some Xq] and by diffeomorphism invari- 
ance, we may assume that xq = 0. Thus in the parabolic region to be characterized 
below, it is natural to study the system in the renormalized variables r and a 
defined by 

(3.1) r(t):=-log(T-*), 
and 

(3.2) a{x7 t):=^L=e^s(x,t), 

where s(x,t) denotes arclength from the point x — 0. Note that a is a natural 
choice for the parabolic region, in which time scales like distance squared. Use of r 
is not necessary but is convenient for the calculations that follow. (Compare [3].) 



The evolution equation satisfied by ip, equation (2.6b I, implies that at the local 
maximum ("bump") s(t) closest to the right pole, one has ipt(s,t) < (n — l)/ip. 
On the other hand, applying a version of the maximum principle to the same 
equation shows that the radius of the smallest neck satisfies the condition V'min(i) > 
y/(n— 1)(T — t). (For the proof, see Lemma 6.1 of [2].) Hence, we introduce the 
rescaled radius 

(3.3) U(a,r):= , ^ . 

y/2(n-l)(T-t) 

Note that U is a function of the coordinate a. In Appendix [Aj we compute the 
evolution equation satisfied by U, which is 



(3.4) U T = U a(T - (| + rdj U a + (n- 

where 

^da. 

o U 

Note that equation (|3.4|) is parabolic, but contains a nonlocal term. 



"*As seen in Section [3] some though not all of the solutions we construct here are reflection 
symmetric. 
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Motivated by Perelman's work |20j in dimension n + 1 = 3, one expects the 
solution to be approximately cylindrical a controlled distance away from the pole, 
depending only on the curvature there. (Essentially, this is contained in Corol- 
lary 11.8 of [20) ; for more details and a refinement, see Sections 47-48 and in 
particular Corollary 48.1 of [19].) This expectation effectively determines the be- 
havior we anticipate in the parabolic region f2 par . 

Ansatz Condition 1. As r /* oo (equivalently, as t /*T), the rescaled radius U 
converges uniformly in regions \a\ < const, where 

U(a,r) 1. 

Since an exact cylinder solution is given by U{ct,t) = 1, the parabolic region 
57 par may be characterized as that within which the quantity 

(3.5) V(a,r):=U(a,r)-l 
is small. 

By substituting U = 1 + V in equation ( |3.4[ )p|one finds that V evolves by 

(3.6) V r = AV + N(V), 
where A is the linear operator 

(3.7) AV:=V trir -~V + V, 
and TV consists of nonlinear terms, 

with / = f° V a<J / (1 + V) da its sole non-local term. 

The operator A appears in many blow-up problems for semilinear parabolic 
equations of the type u t = u xx +u p (e.g. [15]), and also in the analysis of mean cur- 
vature flow neckpinches (e.g. [4]). The operator A is self adjoint in L 2 (R, e~ a / 4 da). 
It has pure point spectrum with eigenvalues {Afc}^L , where 

(3.9) Xk: = l -\- 

The associated eigenfunctions are the Hermite polynomials hk(a). We adopt the 
normalization that the highest-order term has coefficient 1, so that ho(a) = 1, 
hi (a) — a, h 2 (a) = a 2 — 2, h 3 (a) = a 3 — 6a, and in general, 

k 

h k(o-) =J2 r U a ^ 

j=o 

for certain determined constants r)j, with r/k = 1. Recall that this sum for the 
Hermite polynomial hk contains only odd powers of a if k is odd, and even powers 
if k is even. 

In the parabolic region f2 par , which may be regarded as the space-time region 
where one expects V « 0, the linear term AV should dominate the nonlinear term 
N(V), as long as V is orthogonal to the kernel of A. (We argue below that this 
orthogonality follows from the Ansatz adopted in this paper.) Therefore, we may 

^See equation 1A.71. 
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assume that V ~ V , where V is an exact solution of the approximate (linearized) 
equation V T = AV. 

The general solution to the linear pde V t = AV can be written in the form 
of a Fourier series by expanding in the Hermite eigenfunctions, namely V(<t,t) = 
Sfclo °ke XkT h k (°~) , where all coefficients b k except possibly b 2 are constants. If we 
treat V as the first term in an approximate expansion for solutions of the nonlinear 



system (3.6|, then we are led to write 



(3.10) V(<t,t) « b e T + b x e T ' 2 h x (a) + b 2 (T)h 2 (<r) + ]T b k e x » T h k (a) , 

fe=3 

where b 2 (T) is allowed to be a function of time, since it corresponds to a null 
eigenvalue of the linearized operator and hence to motion on a center manifold 
where the contributions of N(V) cannot be ignored^] In a standard initial value 
problem, the coefficients b k are determined by the initial data. Here, they may be 
determined by matching at the intersection with the adjacent regions. 
We now make a pair of additional assumptions. 

Ansatz Condition 2. (i) The solution does not approach a cylinder too quickly, in 
the sense that sup^^gQ |V(<t,t)| > e~ Cr for some C < 00. (ii) The singularity 
time T < 00 depends continuously on the initial data go- 

The key consequence of Ansatz Condition [2] which we explain below, is that 
one eigenmode dominates in the sense that 

(3.11) V(*,t) « V k (a,r) := b k e XkT h k {a) . 

Because A has pure point spectrum in L 2 (M, e~ a2 / 4 da), this is reasonable and 
consistent with Ansatz Condition [T] and with the fact, proved in [2], that singu- 
larity formation at the pole is nongeneric in the class of all rotationally symmetric 
solutions. 



In fact, we may assume that equation (3.11 ) holds for some k > 3. Here is why. 
Part (i) of Ansatz Condition [2] implies that some b k ^ 0. (Compare Section 2.16 
of [3].) Part (ii) implies that one can without loss of generality choose the free 
parameter T = T(go) so that &o — 0. (This too is true for the nondegenerate 
neckpinch; see Section 2.15 of [3].) Ansatz Condition [l] implies that b\ — 0. In 
Section [6j we discard k = 2 as a consequence of Ansatz Condition |1J Hence we 
have k > 3, which implies orthogonality to the kernel of A, as promised above. 

We can now obtain a rough estimate for the size of the parabolic region. Re- 
call that the main assumption made for the parabolic region, Ansatz Condition [T] 
applies only so long as U rj 1, hence as long as 

(3.12) V{a, t) = b k e XkT a k + 0(e X " T cr k - 2 ) (a -> 00) 

is small, hence only as long as \a\ -C e^~^ T . We label the region characterized by 
larger values of |er| the intermediate region. 

4. The intermediate region 

The assumptions we have made to obtain a formal expansion in the parabolic 
region are expected to be valid only so long as V(a, r) is much smaller than one. 

^Recall that 62(7") = (8t) — 1 for the rotationally symmetric neckpinch [3]. 
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We thus leave the parabolic region and enter the next region for those values of 
(<r, r) for which V(cr, r) is of order one. 

By Ansatz Condition [2] and its consequence (3.11), the dominant term in the 
Fourier expansion of V is a polynomial with leading term comparable to e XkT a k . 
So we find it useful to introduce a new space variable 

(4.1) p:={e XkT ) 1 ' k a = e^-^ T a, 

and use p to demarcate the intermediate (p « 1) region. 
In the intermediate region fi; nt , we define 

(4.2) W(p,t) = U(a,r) ( = 1 + V(a,r)). 

Recalling the evolution equation (3.4) satisfied by U, one computes that W satisfies 
the PDE 

\{W-W~ x )-^W p 

= e-W t - ^ + ( „ _ ^ _ n J> dp | . 

This appears to be a small time-dependent perturbation of an ODE in p. We make 
this intuition precise in the following additional assumption: 

Ansatz Condition 3. Wt, W p , and W pp are all bounded in the intermediate region. 

Ansatz Condition M implies that %W P — \ W + \ — o(l) as r -» oo. This 
suggests that W(p,t) w W(p,t), where W solves the first-order ODE 

(4.3) ^w p -^W+^W- x =0. 



One readily determines that the general solution to equation (4.3) is 

(4.4) W(p,t) = y/l-(p/c)". 

Notice that in solving this ODE, we obtain a "constant of integration" c(t) which 
may a priori depend on time. Matching considerations will show that c is in fact 
constant. (The minus sign above is chosen so that the singularity occurs at the 
right pole if c > 0.) 

In order to match the fields in the intermediate region with those in the par- 
abolic region, it is useful to consider the asymptotic expansion of W about p = 0; 
one gets 

W{p,t) = l-\(p/c{t)) k - l -{p/c(t)Y k ---- . 
Thus for < p = (e XkT ) 1/k a < 1, Ansatz Condition^ implies that 

(4.5) W(p, t) » W(p, t) = 1 - ^e A * T a k + • • • , 



which matches the parabolic expansion (3.12) if c(t) is determined by bk, which is 
constant for all k > 3. It is convenient in what follows to regard bk as determined 
by c > 0; thus one has 

(4.6) bk = -\c- k . 
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Going in the other direction, we expect to transition to the tip region where 
W ~ 0, namely where p ~ c. It is unsurprising that one gets the same conclusion 



from equation (3.12| by solving V ~ —1 for large \a\. 

5. The outer region 



For p » 1, which marks the outer region Q ou ti we unwrap definitions (3.2 1, 



(3.3), (4.1), and (4.2) to obtain 



(5.1) V(s, t) 2 = 2(n - l)e- T W(e T/k s, r) 2 « 2(n - 1) 



<*->-(;)' 



If k > 3 is odd, then equation (5.1) implies that the <?(i)-measure of the open 
set (0, 1) x §" goes to zero as t /*■ T. Moreover, the t = T limit profile of i\) can 
vanish to arbitrarily high order, behaving like s k l 2 as s /* 0. 



If k > 4 is even, then equation (5.1) implies that the solution encounters a 



global singularity at t — T, with the whole manifold shrinking into a non-round 
point. 

In either case, a set of nonzero <7(0)-measure is destroyed, in contrast to the 
nondegenerate Type-I neckpinches studied in [2,3, which become singular only on 
the hypersurface {xq} x Note that it is possible for local Type-I singularities to 
destroy sets of nonzero measure. See examples by one of the authors |12) . as well 
as recent work of Enders, Miiller, and Topping [11 . 

6. The tip region 

Obtaining precise asymptotics at the right pole is complicated by the fact that 



the natural geodesic polar coordinate system (2.1) becomes singular there. As in 
[Tj. we overcome this difficulty by choosing new local coordinates. 

The fact that ^ s (s(l, t), t) = —1 for all t < T implies that ip s < in a small 
time-dependent neighborhood of the right pole x = 1. So we may regard i>{s, t) as a 
new local radial coordinate, thereby regarding s as a function of ip and t. (Compare 
0] and pp.) More precisely, there is a function y(tp,t) < defined for small ip > 
and times t near T such that 

(6.1) ^ s {s,t)=y{^{s,t),t). 
In terms of this coordinate, the metric takes the form 

(6.2) g = y(4>{ S , t), ty 2 (d^) 2 + ^ 2 3can , 

with y(i/),t) the unknown function whose evolution controls the geometry near the 
pole. 

For convenience, we replace y < by the quantity z = y 2 , which evolves by the 
pde z t = ^[z], where 

(6.3) := \t\) 2 zz^ - ^(■ipz^,) 2 + (n - 1 - z)ipzy + 2(n - 1)(1 - z)z j . 

For the tip region, we now introduce a i-dependent expansion factor for the 
radial coordinate, setting 

(6.4) 7 ( S ,t) : =r(T(t))V(M), 

with r to be determined below by matching considerations. Defining 

(6.5) Z(7,t):=*fcM) ( = y 2 0M)), 
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one determines from equation (6.3) that Z satisfies 

(6.6) r- 2 (z t + r- 1 r T7 z 7 ) = y 7 [z], 



where is the operator appearing in equation (6.3), with all ip derivatives re- 
placed by 7 derivatives. 

As we observe in Appendix [Sj the ode £F 7 [Z] = admits a one-parameter 
family of complete solutions satisfying the boundary conditions Z(0) — 1 and 
Z(oo) = 0. Any such solution is given by 



where a > is an arbitrary scaling parameter, and B is the profile function of the 
Bryant steady soliton metric, 

9B = B(r)- 1 (dr) 2 + r 2 g can , 

discovered in unpublished work of Bryant, who proved that it is, up to homothety, 
the unique complete rotationally symmetric non-flat steady gradient soliton on 
]R n+1 for n > 2. Recent results of Cao and Chen (followed by a simplified proof by 
Bryant) allow one to replace the words "rotationally symmetric" in the uniqueness 
statement above by the words "locally conformally flat" 6 . Uniqueness under the 
assumption of local conformal flatness for n + 1 > 4 was also proved independently 
by Catino and Mantegazza. Recent progress by Brendle allows one to replace local 
conformal flatness by the condition that a certain vector field V := Vi? + g(R)Vf 
decays fast enough at spatial infinity [5]. (Here, / is the soliton potential function, 
and g(R) is chosen so that V vanishes on the Bryant soliton.) We briefly review 
some relevant properties of the Bryant soliton in Appendix |B| 

Numerical studies for Ricci flow of rotationally symmetric neckpinches [131 114] 
strongly support the contention that for certain initial data, the flow near the tip 
approaches a Bryant soliton model. Results of Gu and Zhu also support this expec- 
tation [16] . Therefore, in the tip region, we adopt the following assumptions, which 



suggest that a solution of (6.6) should be approximated by a steady-state solution 
Z of JylZ] — 0, that is, by a suitably scaled Bryant soliton. These assumptions 
complete our general Ansatz. 



Ansatz Condition 4. In the tip region, we assume that the lhs of equation (6.6) 
is negligible compared to the RHS fort m T. To wit, we assume that (i) T(r) ^> e T I' 1 , 
(it) T r = o(T 3 ), and (in) Z t = o{T 2 ), all as t /• T. 

As noted above, the choice of the expansion factor T(t) is determined by match- 
ing at the intersection of the tip and parabolic regions. We now discuss this de- 



termination. Below, we verify that the form of T(r) we obtain in equation ( 6.11 ) 
satisfies the conditions of Ansatz Condition [2] as long as k > 3. 

To study the consequences of matching at the tip-parabolic intersection, we 
first recall that equation (3.11), with a = e T l 2 s, implies that 



(6.7) V(M) = v / 2(«-l)e ( - 1/2+Afc)r [l + o(l)]6 fe a' £ 
and 

(6.8) ip„(s, t) = v / 2(n-l)e XkT [l + o(l)]fc6 fc cr fc - 1 

hold for large \a\. On the other side of the interface — at the outer boundary of the 
tip region — it follows from Proposition [l] in Appendix |B| and the choice (c 2 = 1) 
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of normalization there that 

(6.9) ^(«,f)«yi(^« or- V" 1 (*>*)■ 



Comparing (6.7), (6.8 1 , and (6.9| at p = c, i.e. at a — e k T c, and using (4.6), 

a 

(6.10) 



one finds that the asymptotic expansions match provided that 

k(n - 1) 



2c 
and 

(6.11) r = e( 1 -^) T . 

This yields a one-parameter family of formal solutions indexed by the scaling pa- 
rameter c > 0. More precisely, the free parameters a, and c obtained in the tip, 
parabolic, and intermediate regions, respectively, are reduced by the tip-parabolic 
and parabolic-intermediate matching conditions to a single parameter, which we 
take to be c. 



The choice of F in equation (6.11) satisfies properties (i) and (ii) of Ansatz 
Condition[4]provided that k > 2. Moreover, the first-order term in our approximate 
solution is stationary, so that property (iii) is satisfied automatically. 

We now recall from Section [2] that an SO(n+ l)-invariant metric has two dis- 
tinguished sectional curvatures, which we call K and L. In terms of the present set 
of local coordinates, the sectional curvatures of g = z^ 1 (dip) 2 + ip 2 g can are given 
by K — —z. l p/{2tp) and L = (1 — z)/ip 2 . At the pole x = 1, these are equal and 
easily computed. Again by Proposition [T] in Appendix [5] one has 

(6.12, A-|„, . L U ,, . to i^|<M . J . jf -fg^ E , 

where k > 3. It follows that these formal solutions exhibit the characteristic Type-II 
curvature blow-up behavior expected of degenerate neckpinch solutions. Notice that 
the limiting behavior as k oo matches the (T — t)~ 2 blow-up rate first observed 
by Daskalopoulos and Hamilton in their rigorous treatment [9] of complete Type-II 
Ricci flow singularities on M 2 , where Ricci flow coincides with the logarithmic fast 
diffusion equation u t — A log?/. The asymptotic profile of its blow-up was derived 
formally by King |18j and recovered rigorously by Daskalopoulos and Sesum [10 . 



7. Conclusions 

Gu and Zhu |16j prove that Type-II Ricci flow singularities develop from non- 
generic rotationally symmetric Riemannian metrics on S" +1 (n > 2), having the 
form g — (ds) 2 + ip 2 {s) <? can in local coordinates on §™ +1 \{P±}. 

Our work above describes and provides plausibility arguments for a detailed 
asymptotic profile and rate of curvature blow-up that we predict some (though 
not necessarily all) such solutions should exhibit. We summarize our prediction as 
follows. 

Conjecture. For every n > 2, every k > 3, and every c > 0, there exist Ricci 
flow solutions g(t) that satisfy the conditions outlined in our Basic Assumption and 
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l:(T-f) 




HT-tf 



Figure 1. An asymmetric degenerate neckpinch viewed at various scales. 



develop a degenerate neckpinch singularity at the right pole at some T < oo. The 
singularity is Type-II — slowly forming — with 

J s T + J Rm(M)l ~ {T-%-m 

attained at the pole. Its asymptotic profile is as follows, where s{x,t) represents arc- 
length with respect to g(t) measured from the location of the smallest nondegenerate 
neck. 

Outer Region: As t /T, one has 



V(M) = [l + o(l)]W2(n-l) 



c -<>-(;)' 



holding for —e < s < c(T — t) 1 ^ if k is odd, and for \s\ < c(T — t) 1 '* if k is even. 
Intermediate Region: As t /* T, one has 



y/2(n-l)(T-t) 

on an interval e(T - t) 1/k < s < e _1 (T - i) 1/fc 
Parabolic Region: As t /•T, one has 



[l + o(l)U 1 



{s/c) k 
T-t 



4>{ s ,t) i + Q (i) V(r-t)fc / s \ 

2c k T-t k \^T~t) 



y/2(n-l){T-t) 

on an interval E\JT — t < s < e(T — t) x / k , where hk(-) denotes the k th Hermite 
polynomial, normalized so that its highest-order term has coefficient 1. 

Tip Region: A Bryant soliton forms in a neighborhood of the pole. Specifically, 
with respect to a rescaled local radial coordinate 



(T-ty- i / k 
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near the right pole, in which the metric takes the form 

g = Z( 1 ,t)- 1 (dyj) 2 + yj 2 g can , 

one has 

z (7 , *) = [i + od)] b(^) asn 

where B denotes the Bryant soliton — up to scaling, the unique complete locally 
conformally fla^non-flat steady gradient soliton on M n+1 . 




Figure 2. An asymmetric degenerate neckpinch viewed without 
rescaling in the (non-geometric) x coordinate system. 

Many aspects of this work are familiar. Indeed, our predicted rate of curvature 
blowup matches that of the examples of Type-II mean curvature flow singularities 
rigorously constructed by Velazquez and one of the authors [4] . The "global singu- 
larities" encountered by the symmetric (k even) profiles considered here agree with 
the intuition obtained from such rigorous examples for mean curvature flow. More- 
over, the Ricci flow singularities numerically simulated by Garfinklc and another of 
the authors [131 114] are qualitatively similar to the case k — 4 considered here. 

On the other hand, it was perhaps not obvious a priori that a Type-II Ricci 
flow solution would vanish on an open set (0, 1) x S" of the original manifold. This 
occurs for the asymmetric (k odd) profiles considered here, which correspond to 
the "intuitive solutions" predicted and sketched by Hamilton [171 Section 3]. 

Motivated by the fact that conjectures provide direction and structure for the 
development of rigorous new mathematics, it is our hope that the formal derivations 
in this paper facilitate further study of Type-II (degenerate) Ricci flow singularity 



Uniqueness also holds if local conformal flatness is replaced by a suitable condition at spatial 
infinity; see [5]. 
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formation. In particular, we intend in forthcoming work to provide a rigorous proof 
that there exist solutions exhibiting the asymptotic behavior formally described 
here. Further study is also needed to determine what (if any) other asymptotic 
behaviors are possible. 

Appendix A. Evolution equations in rescaled coordinates 

In this appendix, we (partially) explain our choices of space and time dilation 
in the parabolic region, and we derive the evolution equations satisfied by a rescaled 
solution. 

Let T < oo denote the singularity time; let 
(A.l) r = -log(T-t); 

and let 

(A.2) a = e 0T s, 

where j3 is a constant to be chosen. Define a rescaled solution U (a, r) by 
(A.3) t) = ^2(n - l)e~ aT U(a, r), 

where a is another constant to be determined. 
It is straightforward to calculate that 

1> t = ^2(n-l)e (1 -^ T (U T + a T U„ - aU), 

i/j s = y/2(n-l)e^-^ T U a , 

As = V2(«-l)e (2 ^ Q)r C/ ffff . 

The factor a T is necessary for a and r to be commuting variables. It is given by 

(A.4) a T = l3a + e - r/2 ^ = (5<t + nl, 

where I is the non-local term 

(A.5) ' 



It follows from equation (2.6b I that U satisfies 

(A.6) eV- 2 n T {U r + a T U a - aU) = XJ aa + (n - 1)^ - e 2 ^ )T ~. 

If f} < 1/2, then for r 3> 0, U should approximate a translating solution of a first- 
order equation. If f3 > 1/2, then for r ^> 0, U should approximate a stationary 
solution of an elliptic equation. The choice j3 = 1/2 is thus necessary if one expects 
U to be modeled by the solution of a parabolic equation^] 

Now suppose that j3 = 1/2 and write U = 1 + V. It follows from the consider- 
ations above that V evolves by 

(O \ V 2 e (2a-l)r 

(A.7) V T = V arT - (_ + nl) V a + (n- 1)^ + a(l + V) - 2(TTVy . 



Caveat: it is not necessarily the case that a parabolic equation will dominate in the "para- 
bolic region;" for example, see the degenerate singularity considered in pQ. 
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Linearizing about V = 0, i.e. U = 1, one finds that V T = AV+ 0<y 2 ), where 



A:V^V a , 



a+-e 



(2q-1)t 



V 



a--e 



(2q-1)t 



The choice a — 1/2 thus results in A becoming the autonomous linear operator 

a , 



(A- 



A : V ^ V aa - -V a + V. 



Appendix B. The Bryant Soliton 

The Bryant soliton, discovered in unpublished work of Robert Bryant, is up 
to homethetic scaling, the unique complete non-flat locally conformally flat steady 
gradient soliton on K n+1 for n > 2. (See [6] and [7].) Uniqueness also holds under 
the assumption that a vector field V := V R+ g(i?)V / decays fast enough at spatial 
infinity. (Recall that / is the soliton potential function, and g(R) is chosen so that 
V vanishes on the Bryant soliton; see [5].) 

Numerical simulations by Garfinkle and one of the authors suggest that a de- 
generate neckpinch solution should converge to the Bryant soliton after rescaling 
near the north pole [131 114) . Results of Gu and Zhu also support this expectation 



Here we recall some relevant properties of these solutions. It is convenient to 
consider a one-parameter family of Bryant soliton profile functions B(-) depend- 
ing on a positive parameter that encodes the scaling invariance mentioned above. 
For more information about the Bryant soliton (including proofs of the following 
claims), see Appendix C of pQ and Chapter 1, Section 4 of [8]. 



Proposition 1 (Properties of the Bryant soliton profile function). 



(1) The ode 3>[z] = 0, where 

3>Lz] := — < r 2 zz rr (rz r ) 2 + (n — 1 — z)rz r + 2(n — 1)(1 — z)z 

r 2 I 2 

admits a unique one-parameter family of complete solutions satisfying 
Z(0) = 1 and Z(oo) = 0. These are given by 



Z{r) = B 



for g > 0, where B is the Bryant soliton profile function. Each member 
of the one-parameter family of complete smooth metrics given by 

g = Z- 1 (r)(dr) 2 + r 2 g c&a , 

is called a Bryant soliton. 

(2) B{r) is strictly monotone decreasing for all r > 0. 

(3) Near r = 0, B is smooth and has the asymptotic expansion 

„, n ,o n ,o a n(n — 1) ,o K 

B(r) = 1 + b 2 r 2 + —Mr* + \ ' fir 6 + ■■■ , 
n + S (n + 3)(n + 5) 

where b 2 < is arbitrary. 
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(4) Near r = +00 , B is smooth and has the asymptotic expansion 

„/x ? 4-n 9 , (n — 4)(n— 7) o fi 

fir = c 2 r- 2 + -c 2 2 r- 4 + 1 > 4^ + • • • , 

n — 1 (n — lj^ 

where C2 > is arbitrary. 

The arbitrariness of 62 and C2 encodes the scaling invariance of the Bryant 
soliton. In this paper, we fix C2 = 1 in order to make explicit the dependence on 
the scaling parameter g > when writing Z(r) = B(r/g). 
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